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1. Experimental data and analysis methods 

This paper presents inclusive spectra of identified hadrons produced in inelastic p+p interac¬ 
tions at 20, 31, 40, 80, 158 GeV/c. Three analysis methods were used: 

• For spectra the method [1], which is based on the fact that the majority of negatively 
charged particles are 71^ mesons. The contribution of other particles is subtracted using 
calculations based on the EPOS model [2], 

• for 71+, 71^, K^, K and p, p the dE/dx method, which uses information on particle energy 
loss (dE/dx) in the TPC gas to identify particles, 

• for 71+, 71^, K^, K and p, p the tof — dE/dx method, which in addition to the dE/dx in¬ 
formation uses time of flight (tof) measurements. Combined tof — dE/dx analysis provides 
excellent separation of different hadron species close to the mid-rapidity region. 

Results were corrected for trigger and event selection biases, contribution of non-target inter¬ 
actions as well as detector inefficiencies and feed-down from weak decays and secondary interac¬ 
tions. 

2. Results on single-particle spectra 

Spectra of transverse mass mj of negatively and positively charged pions, kaons, protons and 
A hyperons produced in inelastic p+p interactions at 158 GeV/c at mid-rapidity are presented in 
Fig. 1 (left). The corresponding data for central Pb+Pb collisions measured by NA49 [3, 4, 6] are 
shown in the right panel of Fig. 1 (right). These data sets were fitted using the simplified blast wave 
model parameterization [7]: = AimjKi ( )/q ( Pr^inhp ^ ^ transverse flow velocity 

pT was calculated from the equation p = tanh^^pj- 

The value of pr is much higher in central Pb+Pb than in p+p reactions. Transverse mass 
spectra are approximately exponential in p+p interactions. The exponential dependence is modified 
by the transverse flow in central (7%) Pb+Pb collisions. Transverse collective flow in Pb+Pb can 
explain the significant spectra shape difference between p+p and central Pb+Pb reactions. 

Rapidity distributions of 71^ produced in inelastic p+p collisions [1] and energy dependence 
of the width a of the rapidity distribution divided by beam rapidity (ybeam) or by the hydrodynamic 
model prediction (Gls) [8, 9] are presented in Fig. 2. The shape of the rapidity distribution is 
approximately Gaussian, however the best fit was obtained by a sum of two Gauss distributions. 
The width of the 71^ rapidity distribution divided by beam rapidity decreases with the collision 
energy, a / Gls and a /ybeam are smaller in p+p than in Pb+Pb interactions. No significant difference 
of the energy dependence is observed for the width of the 7i^ rapidity distribution in p+p and Pb+Pb 
interactions. 

The results presented in this paper are an important step towards the study of the system size 
dependence of the signals of the onset of deconfinement observed in central Pb+Pb collisions, the 
’’kink”, ’’horn” and ’’step” [4, 5]. The 7i multiplicity at SPS energies shown in Fig. 3 increases 
faster in central Pb+Pb than in inelastic p+p collisions (kink) and the two dependencies cross each 
other at about 40A GeV/c. 
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Figure 1: Transverse mass spectra at mid-rapidity measured by the NA61/SHINE experiment in inelastic 
pH-p interactions (left) and corresponding results of NA49 results from central Pb-tPb collisions at 158A GeV 
(right). 
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Figure 2: rapidity distributions measured in inelastic p-tp collisions at 20, 31, 40, 80 and 158 GeV/c 

(left). Dependence of the scaled width of the rapidity distribution on collision energy (center and right). 
World data are taken from Refs. [10, 11, 12]. Results for the width of rapidity distributions are not corrected 
for isospin effects. 


Preliminary NA61/SHINE measurements allow to significantly improve the world data on the 
inverse slope parameter T of transverse mass spectra of kaons [14] as well as the ratio of and 
71+ yields at mid-rapidity [15, 16]. 

The mj spectra were parametrized by the exponential function: 


d^n Spt 

dprdy T^ + iukT ^ 


\Jp\+ m\ - rriK 


T 



(2.1) 


where S and T are the yield integral and the inverse slope parameter, respectively. The S and T 
parameters were fitted to spectra of A'+ and K mesons at mid-rapidity as shown in Fig. 4. 
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Figure 3: Energy dependence of n multiplicity as a function of collision energy. ALICE points were derived 
from Refs. [11, 13]. 




Figure 4: Transverse momentum spectra of and K mesons produced at y « 0 in inelastic p+p interac¬ 
tions measured by NA61/SHINE. 


Figure 5 presents the energy dependence of the inverse slope parameter T of kaons. Surpris¬ 
ingly the NA61/SHINE results from inelastic p+p collisions exhibit rapid changes like observed 
in central Pb+Pb interactions. World data for p+p and Pb+Pb/Au+Au reactions are plotted for 
comparison and were taken from Refs. [14, 17, 18, 19]. 

The acceptance of the identification methods used at the moment by NA61/SHINE does not 
allow to measure spectra of positive pions close to mid-rapidity. Negatively charged pion spec¬ 
tra were measured by NA61/SHINE [1] using the method. This method cannot be used for 
positively charged pions due to the large contribution of protons and positive kaons. 

Therefore, in order to determine the mid-rapidity yield of positively charged pions the follow¬ 
ing procedure was used: 

• The ratio of measured 71+ and yields was calculated (within the acceptance of the tof — 
dE/dx and dE/dx methods) and compared with model predictions as a function of collision 
energy, see Eig. 6. The agreement between the data and the EPOS model [2] predictions is 
better than 0.1%. 
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Figure 5: Energy dependence of the inverse slope parameter T of transverse mass spectra of kaons for 
inelastic p+p interactions and central Pb+Pb/Au+Au collisions. 


The mid-rapidity ratio of the 71+ and tt yields was calculated from the EPOS model. 

The 71+ mid-rapidity yield was calculated as the product of the measured 7i^ yield at y = 0, 
the measured ^ ratio in the tof — dE jdx acceptance and the EPOS correction factor: 


71 ^ 


71+ (y = 0) = TT (y = 0) — {tof - dE/dx) Cmc, 




( 2 . 2 ) 


where: 


Cmc — 


i^(y = o) 

^ {tof) 


:5%. 


(2.3) 


J MC 


The final correction factor as a function of energy is shown in Pig. 6. 


The energy dependence of the A'+/7r+ and ratios at mid-rapidity for inelastic p+p in¬ 

teractions and central Pb+Pb/Au+Au collisions is presented in Pig. 7. Surprisingly the NA61/SHINE 
data suggest that even in inelastic p+p interactions the energy dependence of the A'+/7r+ ratio 
exhibits rapid changes in the SPS energy range. However, the horn structure is significantly re¬ 
duced/modified in comparison to that observed in central Pb+Pb collisions. Additionally world 
data [15, 16, 20, 19, 13] is plotted to establish the trend outside the SPS energy range. The 
NA61/SHINE results were compared to theoretical models in Ref. [25] , namely EPOS [2], UrQMD [22, 
21], Pythia 8 [24] and HSD [23]. Pigure 8 shows that these models do not describe well the 
NA61/SHINE results on p+p interactions. The high precision of the measurements should allow 
for significant improvement of the models. 

Proton transverse momentum spectra and mean transverse mass {mj) around mid-rapidity in 
inelastic p+p interactions at SPS energies are presented and compared with theoretical models in 
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Figure 6: The ratio of 7r+ and yields in inelastic p+p interactions as a function of collision energy. 
The ratio was calculated in the acceptance of the tof — dE/dx and dE/dx identihcation methods. The 
NA61/SHINE results are compared to predictions of the EPOS model (left). Correction factor used to obtain 
;r+ yields from the measured yields (right), see Eq. 2.3. 
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Figure 7: The /tC and K jn ratios in inelastic p+p and central Pb+Pb/Au+Au reactions plotted versus 
collision energy. 


Fig. 9. Mean mj of protons was calculated from the data parametrization: 


S PT 


d^n 


dprdy T^ + tUpT 


exp 


fpT+mj-mp 


(2.4) 


It increases slowly with collision energy. Neither UrQMD nor HSD describe this behavior [25]. 


3. Conclusions 

The NA61/SHINE measurements show that p+p collisions are unexpectedly interesting. Sur¬ 
prisingly the NA61/SHINE results suggest that even in inelastic p+p interactions the energy depen- 
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Figure 8: Comparison of K/n ratios in inelastic p+p collisions at SPS energies with theoretical models. 
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Figure 9: Proton transverse momentum spectra (left) and mean transverse mass mj (right) around mid¬ 
rapidity in inelastic p+p interactions at SPS energies compared with theoretical models. 

dence of the A'+/71+ ratio and the inverse slope parameter of kaon transverse mass spectra exhibits 
rapid changes in the SPS energy range. However, the structures, e.g. the step and horn, are signif¬ 
icantly reduced/modified compared to those observed in central Pb+Pb collisions. High precision 
NA61/SHINE data present challenges for models and should allow their improvement. 
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